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Введение. Актуальность выбранной темы исследования обусловлена необходимостью проведения быстрой 
оценки состояния и надежности материалов, применяемых в различных конструкциях. Целью работы явилось 
исследование параметров влияния дефекта на отклик поверхности среды на ударное воздействие. ешение 
обратной задачи о восстановлении радиуса дефекта базируется на основе сочетания расчетного подхода и 
применения искусственных нейронных сетей (ИНС). Авторами построена методика восстановления параметров 
дефекта на основе применения численного моделирования и ИНС. 

Материалы и методы. Поставленная задача решается в плоской постановке с помощью метода конечных 
элементов (МКЭ). В работе использовались линейные уравнения теории упругости с учетом диссипации 
энергии. В качестве метода решения краевой задачи применялся метод конечных элементов, реализованный в 
пакете 4А/№5У5. В качестве моделирования процесса применения (ИНС) — комплекс МАТЕАВ. 

Результаты исследования. Разработана конечно-элементная модель слоистой конструкции в плоской 
постановке задачи в пакете А№57У5. Решена задача определения нестационарных колебаний при импульсном 
нагружении для различных вариаций радиуса дефекта. Применено позиционное сканирование объекта 
исследований. Построены графические зависимости амплитуд колебаний точек на поверхности от радиуса 
дефекта. 

Обсуждение и заключения. В качестве результатов исследования зависимостей откликов колебаний от радиуса 
дефекта был разработан подход к восстановлению данного параметра в плоской конструкции на основе 
сочетания метода конечных элементов и ИНС. Исследование показало, что использованного объема данных 
достаточно для успешного обучения построенной модели ИНС и проведения идентификации скрытого дефекта 


в конструкции. 
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Введение. Существует обширный перечень различных строительных, мостовых конструкций, зданий, 
оснований, грунтов, композитных материалов, которые эксплуатируются в различных условиях. Конструкции 
могут иметь достаточно большой ресурс эксплуатации и дефекты различной конфигурации: трещины, полости, 
включения, коррозионные изменения структуры, вмятины и пр. Вследствие этого могут возникать 
непредвиденные предаварийные состояния и отказы конструкции. Для предупреждения аварийного состояния 
требуется решение задачи их ранней диагностики. 

Этот факт определяет необходимость использования различных методов неразрушающего контроля [1— 
2]. Некоторые направления развития отдельных методов диагностики наличия дефектов в различных объектах 
представлены в [3-10]. 

Основа методов состоит в сборе и анализе определенных параметров конструкции, являющихся 
откликом при ее нагружении. В качестве одного из таких параметров может быть вибросмещение 
определенных точек поверхности при импульсном нагружении колебаний конструкции. При этом от места 
возбуждения колебаний распространяется волновой фронт, возбуждающий смещение точек и колебания всех 
элементов конструкции. Задача может быть сведена к рассмотрению измерения скоростных параметров 
распространения волнового поля. Основываясь на анализе публикаций в иностранной литературе, можно 
констатировать, что экспериментальные исследования показывают явную тенденцию к оперированию методом 
«измерения скорости» [11-12]. 

Применение специализированных приборов, которые располагаются в определенных точках объекта и 
осуществляют сбор информации о различных параметрах колебаний, описаны в [13, 14]. Данная методика 
предполагает позиционное, временное и частотное сканирование объектов исследования. Наиболее 
привлекательным является использование акустических датчиков и приемников, установленных на внешней 
поверхности конструкции. Данные датчики могут фиксировать отраженные сигналы от дефектов. Применение 
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диагностических систем и комплексов требует разработки программного обеспечения (1/0), которое позволяет 
производить первичный анализ сигнала. В случае «тонкой» приборной настройки, например, при применении 
специальных алгоритмов обработки по анализу отраженного сигнала, ГО позволяет идентифицировать 
дефектную область конструкции. Соответствующее ГПО может быть разработано на основе использования 
искусственной нейронной сети (ИНС) [15]. Применение ИНС в задачах реконструкции поврежденного 
состояния элементов конструкций описано в работах [16—26]. Применение различных архитектур и алгоритмов 
ИНС описано в [16-20]. Определению дефектов в анизотропных пластинах с помощью ИНС посвящена работа 
[21]. В работе [22] авторы указали преимущества методов идентификации, не требующих предварительного 
построения математической модели объекта исследования. 

В настоящей работе разрабатывается метод реконструкции диаметра круглого дефекта, 
расположенного в полуплоскости слоистой конструкции. Такого рода дефекты могут часто встречаться при 
вскрытии верхних слоев грунтов в виде уплотнений каменных структур под слоями дорожной одежды. Модуль 
упругости и плотность таких структур могут в разы отличаться от основных параметров слоя. В 
математической постановке проблема сводится к решению обратной геометрической задачи теории упругости 
[27]. Рассматривается построение алгоритма восстановления информации об объекте, локализованном в 
определенном месте. Для этого рассматривается применение метода КЭ, анализ результатов поверхностных 
волн и корреляционный анализ зависимости параметров дефекта от волнового поля отклика с применением 
ИНС. 

Материалы и методы. Предлагается подход целью которого является восстановление информации о 


параметрах дефектной области конструкции на основе анализа отклика поверхности среды на ударное 
воздействие. Анализ осуществляется на основе построенных алгоритмов при сочетании расчетного подхода и 
применения ИНС. В качестве образца рассмотрен слоистый материал. 


Нагружение конструкции импульсным ударным воздействием производится на некотором расстоянии 
от исследуемой зоны, а датчики, регистрирующие колебания, располагаются в определенных позициях и 
последовательности в линейном направлении зоны исследования. Могут измеряться поперечные и продольные 
колебания. Интерес представляет период времени, когда волны, отраженные от концов моделируемого участка 
конструкции, не успевают достичь датчика. Таким способом имитируется моделирование реальных условий 
эксплуатации слоистой конструкции грунтового основания. Анализ измеренных амплитудно-временных 
характеристик (АВХ) показывает возможность их использования в обратных задачах восстановления дефектов. 

В качестве инструмента решения обратной задачи реконструкции параметров дефекта используется 
ИНС, которые изначально проектировались для решения задач определения нелинейных зависимостей в 
многомерных массивах данных. В отличие от других алгоритмов, ИНС не программируются, а обучаются на 
множестве данных для исследуемых параметров конструкции. Обучающие выборки строятся путем решения 
прямых задач в конечно-элементном (КЭ) пакете АМ5У5. Обученная сеть, получив новые, неизвестные ранее 
результаты анализа, способна корректно распознать параметры дефекта. 

Эффективные применения аналитического и численного моделирования, достаточно хорошо 


коррелирующие с восстанавливаемыми параметрами элементов конструкций, представлены в [28-34]. 
Постановка задачи. Поставленная задача решается в плоской постановке с помощью метода конечных 
элементов (МКЭ). В настоящей работе используются линейные уравнения теории упругости с учетом 


диссипации энергии, принятой в пакете 4А№5 75 [34, 35]. 


Для упругой среды: 
ри, + ари, о, =Д; 
о; = са (Ев + ВЕ); 
р Ик т ы 
2 

где р — плотность материала; и; — компоненты вектора-функции перемещений; 0;; — компоненты тензора 
механических напряжений; }. — компоненты вектора плотности массовых сил; =; — компоненты тензора 
деформаций; с; — компоненты тензора четвертого ранга упругих модулей; а — неотрицательные 


коэффициенты демпфирования (в АМЗУ5). 
Рассмотрим подробнее механические краевые условия. При определении механических граничных 
условий граница тела представляется в виде объединения непересекающихся областей 5 = 5,5+.5 и» на которых 


задаются следующие условия: 
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— условие закрепления границы или заданных перемещений: 
— 0 
и: |5, =; , 
— условие силового воздействия, при которых заданы компоненты вектора механических напряжений: 
[м — отл, И 
—щ условие гладкого контакта с абсолютно твердым телом — равенство нулю нормальных перемещений и 


Описание модели. Рассматривается слоистая конструкция. Верхний слой жестко сцеплен с 
подстилающим полупространством (рис. 1). В слое |1 расположен дефект в виде круглой конфигурации с 
центром на определенной глубине. Свойства материала представлены в таблице 1. Дефект расположен на 
глубине 7У,,. = 1,5 м от поверхности и на расстоянии Х». = 2 м от точки Р, приложения единичного 
импульсного Р, (рис. 1) нагружения. Импульсное нагружение зависит линейно от времени приложения 
нагрузки (т = 0,003 с), что соответствует параметрам реального ударного нагружения (рис. 2). В качестве 
входного параметра рассматривается поле смещений поверхности в результате кратковременного импульсного 
воздействия. Радиус дефекта варьируется в расчетах как: А; = 0; 0,25; 0,3; 0,35; 0,4; 0,45; 0,5 м. В качестве 
базовой цели ставится определение функциональной зависимости радиуса дефекта от откликов колебаний, 
измеренных в определенных точках конструкции на основе применения ИНС. 


Таблица 1 
Параметры слоев 


Коэффициент 
Наименование Толщина, м Е, Гпа у 
затухания 


2000 00 


| 


х =2 
Дефект К; — 1 2000 0,33 0,001 
У сектЕв= —1,5 


Слой 1 


= Е, Р1. У1, [9 
её 


Ра: Ра: Уа: а 


Подстилающий слой Е\, р1, \1, & 





а) Ь) 


Рис. 1. Модель слоистой конструкции с круглым внутренним дефектом: 
а — описание области с дефектом; 6 — схема модели слоистой конструкции 


РЁ 
р 





Рис. 2. Зависимость импульсной нагрузки Р от времени нагружения 
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КЭ моделирование. Моделирование плоской конструкции (ширина — 10 м) с дефектом проведено в 
КЭ комплексе А№575: для моделирования применялись элементы типа РГАМЕб2 с треугольной и 
четырехугольной конфигурацией с плоской деформацией (3950 узлов и 1890 конечных элементов). Внешний 
слой конструкции (подстилающий) имел большие коэффициенты демпфирования, но сохранял основные 
свойства слоя |. Анализируемое время выбиралось таким, чтобы импульс от крайних границ конструкции не 
возвращался обратно. Это достигалось фиксацией времени резкого возрастания амплитуд колебаний у внешней 
границы слоя 1. Ударная нагрузка прикладывалась в точке Р, на поверхности конструкции. 

Вследствие импульсного нагружения возникают поперечные и продольные волновые смещения точек 
как на поверхности конструкции, так и по всему ее объёму. Пример формы поперечных колебаний конструкции 
в момент времени {Е = 0,01 с представлен на рис. 3. Контрольные точки (40 точек) расположены на поверхности 
на расстоянии АК = 0,1 м друг от друга (рис. 1). Имитация распределения точек на поверхности отображает 
реальную картину измерений. Первая точка измерений расположена на расстоянии 0,1 м от места ударного 


воздействия. 


ШИ 0-128Е-7 
ии 0.106Е-8 
пи -0.106Е-7 
и -0.223Е-7 
ии -0.339Е-7 
(ии -0.456Е-7 
ШЕИ -0-572Е-7 
ии -0.689Е-7 
ШЕИ -0.805Е-7 
ШИ -5.222Е-7 








Рис. 3. Волновые смещения по толщине конструкции при [ = 0,01 с 


На рис. 4, 5 представлены результаты расчетов поперечных смещений (направление ОУ) в различных 
точках поверхности. В процессе решения задачи вычисляли поперечные смещения в контрольных точках 
конструкции на временном интервале # = 0-0,1 с. Особый интерес представляет поле смещений точек 
поверхностного слоя в промежутке времени [ = 0-0,04 с. В данный период времени образуется первичная 
форма прогибов слоистой конструкции и отсутствуют отклики колебаний от возбуждения внешних слоев 
конструкции. 

На рис. 5 приведены максимальные значения амплитуд колебаний в зависимости от времени их 
регистрации. Анализ показывает, что они имеют экспоненциальный приближенный тренд. Для одинаковых 
позиций измерения максимумов амплитуд имеются расхождения амплитуд колебаний в зависимости от 
величины радиуса дефекта. 

На последующем этапе вычислены относительные значения смещения амплитуд ДО, (а), скорости 
АГ, (Ь) и ускорения ДА, (с) для точек №,=1-40. Данные относительные значения вычислялись как разница 
между текущим показателем максимальной амплитуды для точки при К; (при 1=2-—7) и амплитудой параметра 
колебаний при размере дефекта А;=0 (отсутствие дефекта в конструкции): 


АЕ 
ЛД Ин -И 
ЛА=А- т 


Ни рис. 6 приведены трехмерные графики, отображающие зависимости относительных величин 
смещения амплитуд ДЦ, (а), скорости ДИ, (Б) и ускорения ДА, (с) для точек №-1-40 на поверхности 
конструкции и варианта радиуса КЮ; дефекта. 
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Рис. 4. Амплитуды смещений Ц, (а), скорости Г, (Б) и ускорения 4, (с) в зависимости от времени для точек 1-40 на 
поверхности конструкции 
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га Рис. 5. Зависимости максимальных перемещений: 

В а — от времени для точек 1-40 на поверхности конструкции для семи вариантов радиусов дефекта; 
— Ъ — для выделенной области амплитуд 
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АТ, т/5-1 05 


Рис. 6. Относительные значения амплитуд смещения ДИ, (а), скорости ДГ, (5) и ускорения ДА, (с) для точек №,=1-40 на 
поверхности конструкции и варианта радиуса А; дефекта 


Центр дефекта располагался под точкой Ро. Анализ данных графиков показывает, что слева и справа в 
окрестностях дефекта возникают перегибы кривых максимальных амплитуд колебаний относительных 
значений перемещения ДО,, скорости ДГ, и ускорения ДА,. Эта тенденция сохраняется при всех вариациях 
радиуса дефекта К;. Этот признак тем самым может являться показателем локации дефекта. Такие эффекты на 
практике возможно достичь путем позиционного сканирования при регистрации и обработке амплитуд 
колебаний. Кривые относительных показателей перемещения, скорости и ускорения в окрестности дефекта 
достаточно хорошо коррелируют с величиной дефекта. 

Применение нейросетевых технологий в задаче идентификации размеров скрытого дефекта. 
Установим связь между скоростью, ускорением и амплитудами колебаний, распространяющимися в слоистой 
структуре, и радиусом дефекта К;. На основе полученных данных на этапе КЭ моделирования формируется 
обучающий набор, с помощью которого обучается построенная модель ИНС. 

В результате численного решения ряда прямых задач в КЭ программном комплексе 4А№575 были 
получены данные для 40 точек на поверхности слоистой структуры при вариациях радиуса дефекта А; с 
ошибкой при вычислении в пределах 3%. Для каждого набора из 40 значений был установлен соответствующий 
радиус дефекта. Подготовленные обучающие вектора в качестве входных значений содержали относительные 
величины амплитуд колебаний ДИ,(№,„,) скорости ДГ,(М№„) и ускорения ДА,(М№,„;), а в качестве выходных значений 
— радиус дефекта. Всего было проведено 100 численных экспериментов. Таким образом, обучающая выборка 
состоит из 100 векторов для каждого типа исследуемого параметра АЦ, АГ, и ДА,,. Все данные, применяемые в 
обучении нейронной сети, нормализованы и находятся на интервале [0, 1]. 

Для идентификации радиуса дефектного включения применялась полносвязная многослойная модель 
ИНС, смоделированная в комплексе Майа. Модель ИНС содержала 1 слой, состояла из 40 входных нейронов и 
1 выходного нейрона. В качестве входной активационной функции выбиралась 51е11014, а в качестве выходной 
активационной функции устанавливалась линейная зависимость. Сеть обучалась с помощью алгоритма 
обратного распространения ошибки на основе оптимизации Левенберга-Марквардта [37]. 

Функция потерь характеризовалась среднеквадратической ошибкой (МЭЕ). Для обучающей и тестовой 
выборки было построено 100 векторов из показателей относительных величин амплитуд колебаний ЛИ,, 
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скорости ДГ, и ускорения ДА,. Найдена корреляции данных параметров к радиусу дефекта. Установлено, что 8, 
10 и 20 эпох обучения для соответствующих показателей достаточны для достижения необходимого уровня 
качества работы ИНС. На рис. 7 приведена зависимость средне-квадратической ошибки от количества эпох 
обучения модели ИНС. Анализ показывает, что при достижении обучения более 8 эпох для всех относительных 
параметров амплитуд колебаний средне-квадратическая ошибка (М5Р) изменяется незначительно. 


Вез{ УайдаНоп РеМогтапсе 1$ 0.0001885 а{ еросп 7 Вез{ Уа!Чайопт Реогтапсе 15 0.00013677 аЁ еросп 5 


10-2 
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Рис. 7. Изменение средне-квадратической ошибки (МЕ) работы ИНС в процессе её обучения на основе значений 
относительных величин колебаний: а — амплитуд ДЦ,, Ь — скорости ДИ, и с — ускорения ДА, 


Таблица 2 
Оценка точности вычислений 


Используемая ИНС 
АО, 


1670 


Тестирование полученных зависимостей ИНС. Для выборочной оценки восстановления показателя 





радиуса дефекта К; были выполнены наборы из параметров амплитуд колебаний и получены их относительные 
значения. Графическое отображение результатов трех тестовых наборов для относительных величин амплитуд 
колебаний ЛО, (а), скорости ДИ, (Б) и ускорения ЛА, (с) представлены на рис. 8. Данные показатели наборов 
подставляли в соответствующую ИНС и вычисляли значения радиуса дефекта. В таблице 2 приведена 
рассчитанная погрешность восстановленного радиуса К дефекта. Анализ показывает, что наименее точным 
является восстановление радиуса по параметрам ускорения. При этом средняя относительная ошибка 


ЬИр://уезшК-допза.га 


восстановления параметра радиуса К дефекта на основе данных амплитуд колебаний ДИ, и скорости ДИ, не 
превосходит 5%. 
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Рис. 8. Набор рассчитанных относительных величин амплитуд 
колебаний ДЦ, (а), скорости ДГ, (5) и ускорения ДА, (с) 


Обсуждение и заключения. Рассмотрена задача о восстановлении радиуса дефекта на основе 
применения имитационного моделирования колебаний в конечно-элементном комплексе А№575 и ИНС в 
комплексе МАТГАВ. Получены результаты нестационарных колебаний полуплоскости с круглым внутренним 
дефектом. Вычислены поперечные смещения в контрольных точках конструкции при распространении волны 
от импульсного нагружения. В результате был разработан подход к восстановлению радиуса дефекта в 
конструкции на основе сочетания метода конечных элементов и ИНС. Наиболее удачная конфигурация 
архитектуры ИНС, состоящей из | скрытого и | выходного слоев, включала 40 входных нейронов и | выходной 
нейрон. 

Исследование показало, что использованного объема данных достаточно для успешного обучения 
построенной модели и проведения идентификации скрытого дефекта конструкции. При этом относительная 
точность определения радиуса дефекта внутри слоя 1, в случае применения значений амплитуд колебаний 
перемещения и скорости, составляет более 99 %. 

Таким образом построенные ИНС-алгоритмы могут успешно применяться для оценки стратификации 
слоистых конструкций при использовании временного и позиционного сканирования при возбуждении 
колебаний из дальней зоны. 
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